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Abstract 

An analytical model is developed to study the spectra of electromagnetic dissoci- 
ation of two-neutron halo nuclei without precise knowledge about initial and final 
states. Phenomenological three-cluster bound state wave functions, reproducing the 
most relevant features of these nuclei, are used along with no interaction final states. 
The ^He nucleus is considered as a test case, and a good agreement with experimen- 
tal data concerning the shape of the spectrum and the magnitude of the strength 
function is found. 
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1 Introduction 



Coulomb excitation reactions serve as one of the most powerful tools for in- 
vestigating excited states of nuclei. An appealing feature of these reactions 
is the clear understanding of the interaction mechanism. In particular, elec- 
tromagnetic dissociation (EMD) of halo nuclei has revealed an anomalously 
large cross sections due to accumulation of electric dipole (El) strength at low 
energy. The first attempts to interpret the phenomenon of low-energy excita- 
tions addressed a notion of the so-called soft dipole resonance [1-5]. However, 
it is still an open question whether this is a true resonant state, which is being 
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observed through the El excitation. For example, intensive studies of the one- 
neutron halo nucleus ^^Be have shown a very significant enhancement of the 
El strength just above the ^''Be + n threshold where there arc no resonances in 
^^Be. For ^He several different three-body approaches [3,6-9] agree that there 
is no low-lying 1~ resonant state in the continuum. An alternative way to ex- 
plain the accumulation of El strength would be to consider a direct breakup, 
where the final states are contimmm states without a resonant state. In the 
direct breakup mechanism, low-energy dipole excitations occur as a natural 
consequence of the small binding energy or correspondingly of the large dis- 
tance between the charge and the center of mass (CM) of a halo nucleus. In 
this way the El strength function of the one-neutron halo nucleus ^^Be can 
be explained using a simple two-body model with a Yukawa wave function 
(WF) for the initial state and plane waves in the final state [10,11]. However, 
for two-neutron halo nuclei, like ^He, ^^Li, and ^^Be, a three-body picture is 
more appropriate but also more complicated. Microscopic three-cluster cal- 
culations of the El strength function were performed for, e.g., ^He [9,12]. 
However, many difficulties and questions caused by incomplete knowledge of 
the cluster dynamics are still to be resolved. In the present paper we develop 
an alternative approach. We construct phenomenological three-cluster bound 
state WFs of two-neutron halo nuclei that behave correctly at large interclus- 
ter distances, reproduce the nuclear sizes, and incorporate the main features 
of the underlying three-body structure. No interaction three-body WFs are 
used to describe the breakup final states. The model allows an analytic cal- 
culation of the strength functions and can serve as a helpful tool to predict 
the Coulomb disintegration spectra of a variety of two-neutron halo nuclei. 
This can help to conduct new experiments and eventually can provide a better 
understanding of the microscopic structure of halo nuclei. An approach of this 
type has previously been used by Pushkin et al. [13] and the present work can 
be viewed as a development of their model. 



In the present paper the model is formulated. As a test case the El strength 
function of ^He is studied and compared with the experimental data of Au- 
mann et al. [14]. In a subsequent paper the model will be applied to ^^Li and 
^^Be. In Sec. 2 our WFs are described, in Sec. 3 the strength functions are 
obtained and discussed, and in Sec. 4 some conclusions are presented. Details 
of the calculation are given in the Appendix. 



For further reference we list here conventional formulae concerning the EMD. 
We consider the case when an initial ground state is the only bound state in a 
system (note that all known Borromean nuclei possess this property) and all 
possible final states belong to the continuum. In the framework of first order 
perturbation theory the energy spectrum for El Coulomb excitation can be 
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written as 



da(El) _ Nei{E*) 167rMS(El) 
dE ~ he 9 dE ^ ' 

Here E* is the excitation energy, E* = Eq + E, where Eq is the binding 
energy, and E is thus the continuum energy, dB{El) / dE is the dipole strength 
function, and N-ei{E*) [15,16] is the spectrum of virtual photons (where the 
actual number dn of virtual photons equals to N-ei{E*){E*)^^ dE). Since the 
spectrum of virtual photons peaks at low energies, the Coulomb excitation to 
low-lying states is favored as far as there exists a low-energy contribution to 
the dipole strength. 

The El strength function can be written as 
dS(El)_ 1 ^ ^ f-^ .... ...m2 



dE 2Ji + 



^ Mi M=-1,0,1 



dTf \{f\M{El, JM)r S {Ef - E) . (2) 



Here dry is the phase space element for final states, A4(E1) is the dipole 
operator, and |/) are the initial state and the final states in the CM 
subsystem which are normalized as follows 

= (/|/') = 5(r,-r;). (3) 

In the case when discrete quantum numbers enter the labelling of continuum 
final states the 5-function notation adopted above implies inclusion of 5- 
symbols. Similarly, the notation / dry in Eq. (2) may imply the inclusion of 
summation over discrete quantum numbers. 

Since halo nuclei exhibit a large degree of clusterization, low-energy excita- 
tions will mainly affect relative motion between the N clusters. The corre- 
sponding cluster El operator is 




M{El) = W — VeZ,(n-i?cm), (4) 



where fi are the cluster positions, and i?cm is the position of the CM of the 
system. In the case of two-neutron halo nuclei only the core will contribute in 
Eq. (4). 

By summing the strength over all final states one obtains sum rules for the 
process. In particular, the non-energy-weighted cluster sum rule reads in our 
case as 



/ 



~ dS(El)^^_ 3 2/ . 



dE Att 



dE^-Zy{ri), (5) 
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where is the charge of the core, is the distance between the core and the 
CM of the whole system, and the average value is calculated over the ground 
state WF. 



2 Wave functions and transition matrix elements 



In the present analysis we adopt the three-body model description of two- 
neutron halo nuclei. The cluster part of the bound state WF, in the CM 
subsystem, is written as an expansion over hyperspherical harmonics (HH) 
(see e.g. [17]) 



KLSl^ 



^ (f , y) = p-'/' Y: XkIs (p) \^kI m ® Os . (6) 



JM 



Here {x,y} is the set of Jacobi coordinates (A.l), and {^,^25} are the cor- 
responding hyperspherical coordinates. The quantity p — {x^ + is the 
hyperradius, and {fl^} denotes collectively five angles parametrizing a hyper- 
sphere with p — const. We use below that 

dxdy^p^dpdQs- (7) 



The HH, T^KLMi^ (^5)) form an orthonormalized complete set. Harmonics with 
the orbital quantum numbers L, Ml are coupled with spin functions of two 
neutrons 9sms to the total momentum J, M. The other quantum numbers 
labelling the HH are the Jacobi orbital momenta Ix, ly, and the hypermomen- 
tum K. Since the WF (6) should be antisymmetric with respect to the valence 
neutrons it includes only terms with even {l^ + S). More details can be found 
in the Appendix. 

For Borromean systems, having no bound subsystems, which include two neu- 
trons as constituents, the hyperradial functions entering the expansion (6) 
behave asymptotically as (see, e.g., [18]) 

X\{p) ^ Cxex.p{-K,op), as p ^ 00, (8) 

where kq is connected to the binding energy via Eq = (/ifi;o)^/(2m), and m is 
the nucleon mass. Thus, for a nucleus with a small binding energy the WF 
has a long tail which is of importance for peripheral reactions such as EMD. 
One could choose a phenomenological bound state WF using the normalized 
hyperradial function 

X^'Hp) = \^exp(-«op), (9) 

together with the HH from expansion (6) which has a predominant weight. 
The single free parameter, kq, is fixed from the binding energy. Such a model 
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WF, incorporating the K — HH, has been used in Ref. [13] for the -^^Li 
case and led to an analytic expression for the El strength function reproduc- 
ing well the shape of existing experimental data. The total WF behaves as 
p-5/2 for small p but it is still normalizable since the singularity cancels with 
the factor in the volume element, Eq. (7). However, using this model the 
WF is overestimated at small p and thus the ^2ko vahic is smaller than the 
true asymptotic constant. This leads to an underestimation of (p^) and con- 
sequently to an underestimation of the size of the system since, in accordance 
with Eq. (A. 2), 

(i?Ls) = + (A- 2)(i?Ls(core)) 

where A is the mass number, and the last term represents the intrinsic size of 
the core. 

In our model, to cure this feature, we will be add an extra hyperradial term of 
the same exponential form to reproduce simultaneously the true asymptotic 
behavior of the ground state (connected to the binding energy) and its size. 
The corresponding normalized function is 

X^^Hp) = c [exp(-«;op) - exp(-Kip)] , 

2k,ok,i{k,o + Ki) (10) 



where c = 



2 



The parameters kq and ki are fixed using experimental values of binding 
energy and rms radius. With the condition that ki > kq we ensure that the 

second term decays faster than the first, and thus the correct asymptotics is 
preserved. Furthermore, the divergence of the total WF at small p has been 
reduced to The particular form of the function (10) has been chosen to 

be able to perform the calculations analytically. 

We shall retain only one or a few terms in the HH expansion (6). This approx- 
imation is motivated by the rapid increase of the multidimensional centrifugal 
barrier as the hypermoment K increases, see e.g. [17]. The predominant terms 
in the expansion usually correspond to the lowest possible value of K that is 
not suppressed by the Pauli principle. 

In the ^He case the above approximations lead us to the following normalized, 
model WF for the initial bound state ( J'^ = 0+) 
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*(f , y) = ^^-^ [ [v^rl (^^5) + (Qs) 

where Wqq + W2q + W21 = 1. 

Here x^^^ is given by Eq. (10) (or alternatively by Eq. (9)). For the HH retained 
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in Eq. (11) realistic hyperradial functions, in the corresponding expansion, be- 
have rather similarly at large p of interest, so that the model with hyperradial 
functions of the same form is permissible. The a-N interaction is repulsive in 
s states due to the Pauli principle, and it is attractive in p states. Because 
of this the K — contribution to the WF of ^He is small, and the lowest K 
value, which is not Pauli suppressed, is i^' = 2. Prom the K = 2 contributions 
to the WF. given by the second and third terms in Eq. (11), only the r2oo 
term involves the s wave NN attraction and at the same time it includes a 
very large component with the p wave a-N attraction. As a result, the term 
with P^oo should have a predominant weight, and to a first approximation we 
may keep only this single term. One can also retain all three terms entering 
Eq. (11) using stable estimates for the weights of the HH components taken 
from theoretical predictions [17], see Table 1. Similar phenomenological WFs 
may be constructed for other two-neutron halo nuclei. Explicit expressions for 
the HH entering the initial state (11) are given in Eq. (A. 7). 

In Table 1 parameters of our ^He bound state WFs thus obtained are listed. 
The experimental values of the two-neutron separation energy = 0.97 MeV 
and the rms radius Rrms = 2.50 fm are adopted for ^2 and ^3. The i?rms = 1-79 
fm for the function is too small as discussed above. 

In Fig. 1(a) our model hyperradial function (10) is compared with the hyperra- 
dial functions obtained from microscopic three-body calculations of ^He [17]. 
(For other halo nuclei such functions are not known with a sufficient confi- 
dence.) Pig. 1(b) illustrates the ranges of p values important for calculating 
the strength function at energies not far away from its maximum. Typical 
integrands entering the matrix elements (ME), Eq. (14), are plotted. They 
include our hyperradial bound state function and hyperradial components of 
final state WFs. In our case the latter components correspond to = 1 and 
K — 3 (cf. below) from which the K — 1 case is shown. The plot ior K — 3 
looks similar. One can sec that at energies not too far from the peak of the 
strength function, at E' ~ 2 MeV [14], the model hyperradial function is on 
average rather close to the realistic ones at p values of interest. (We note in this 
connection that in the peak region the contribution of the K = component 
of the WF to the net result is much higher than its relative weight.) At ener- 
gies in the maximum region our model ground state WF leads to a strength 
function which is very close to that obtained with the realistic ground state 
WF. Away from the maximum region one should expect more strength with 
our model than with the realistic ground state WF provided that the final 
state WFs are the same. 

We calculate the strength function disregarding the final state interaction 
(FSI). We shall see that a reasonable agreement with experiment both in the 
energy dependence of the strength function and in its magnitude emerges in 
this way. 
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To obtain a differential cross section of EMD one would have to use final state 
WFs with given momenta, including angular information. When the FSl is 
disregarded these WFs are three-body plane waves. To carry out the calcula- 
tions these plane waves could be expanded in products of HH in coordinate 
and momentum spaces, see e.g. [19]. However, in our inclusive case, when we 
arc only interested in the energy dependence of the cross section, we do not 
need directions of the momenta. Thus, instead of using plane waves, we will 
use a set of no interaction final states that include just the coordinate space 
HH. Indeed, the strength function (2) does not depend on the choice of final 
state set (for the same Hamiltonian) provided that the proper orthonormal- 
ity conditions (3) are fulfilled. The states wc choose are labelled by energy, 
hyperspherical quantum numbers L, l^, ly, and quantum numbers J, M, S 
of total momentum and of spin of the valence neutrons. The corresponding 
configuration space WFs include Bessel functions, and they are of the form 



(12) 



Alike the plane waves these functions are solutions to the free-space, six- 
dimensional Schrodinger equation. The continuum energy is related to k via 
Ef = (/?A)^/(2m). see e.g. [19]. The phase space element needed in Eq. (2) to 
integrate over the contributions of states (12) is 

/dr/= E /-'d«: = 4(g)' E J EjdEf (13) 

where the sum is over final states. In accordance with Eqs. (3) and (13) the 
states (12) are normalized to k.~^6{k, — k,') times the ^-symbol with respect to 
the discrete quantum numbers. 

The transition ME between components of the initial and final states with 
given K values obey the selection rule AK = ±1 (the reason is that the 
transition dipole operator is a HH with K — 1). HH have parities (—1)^ so 
that, e.g., in the usual case when the cluster part of the ground state WF 
has positive parity, final states (12) with odd K contribute to the result. 
Taking Eqs. (7), (10) and (11) into account together with the proportionality 
of the transition operator to p, the ME of A^(El) operator (4) will include 
hyperradial integrals of the form 

{k/k,)-'' / J^+2((«M)^)^'/'e-^ dz, (14) 
■/o 

with i = and 1. Note that the whole energy dependence of the total ME is 
included in these integrals. 
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3 El strength functions for two— neutron halo nuclei 

Before presenting more detailed results we shall obtain an approximate general 
property of the strength functions. Let us compare to each other the strength 
functions for a pair of two-neutron halo nuclei. We assume first that the initial 
ground states of both nuclei are dominated by the HH with K = 0. Then the 
largest contributions to the strength functions come from the final states with 
K — 1. At moderate energies predominant contributions to our ME come 
from high p values. For time being let us suggest that the p dependence of the 
ground state WFs at such p values is described by Eqs. (6), (8). As a result 
the energy dependencies of both strength functions will mainly be determined 
by the quantities (14) with i = and K = 1. Thus these energy dependencies 
differ only in the parameter kq. We come to the conclusion that if the strength 
functions of different two-neutron halo nuclei are plotted on the scales E/Eq 
their forms should be similar to each other. This feature is characteristic of 
the no-FSl approximation but, in view of the results below, one can hope 
that it remains approximately valid also beyond this approximation. Another 
limitation is that in reality the high p behaviour (8) may not (in some cases) 
be achieved at p values of interest, and power corrections in the corresponding 
expansion 

CAexp(-/top)(l + ap~^ + . . . ) 

can still be important. Some of these corrections violate the scaling property. 
However, within the ranges of p values effectively contributing to (14), see 
Fig. 1(b), the change of the power terms is relatively small as compared to 
that of the exponential and to a certain approximation the scaling property 
remains valid. 

This property should hold true also when the ground states are dominated by 
a HH with K ^ which is the same for both nuclei, say r2oo with K — 2. In 
such a case the same combinations of hyperradial integrals (14) will define the 
strengths. The scaling property also fulfills in all cases when the contribution 
of the hyperradial integrals (14) with K = 1 dominates the strengths. 

In Fig. 2 several experimentally measured strength functions for ®He and -"^^Li 
are compared in the described way. The spectra are plotted as functions of 
the parameter x = E/Eq. The similarity of these scaled strength functions 
is evident. Such plots may be useful for quick predictions of shapes of El 
strength functions for two-neutron halo nuclei or for checking hypotheses on 
their three-body structure. Shown in this figure is also the analytical strength 
function from a two-body calculation where a Yukawa WF has been used 
together with plane waves in the final state [10,11]. The two-body model 
gives a peak at too low energy thus clearly indicating the importance of using 
a three-body approach. 
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Now let us obtain the El strength function of Eq. (2) in our model. The final 
states (12) which can be reached from the groimd state (11) via the dipole 
excitation are those with K = 1 and K = 3. Other selection rules are Al^ = 0, 
Aly — ±1, and AL = 0, ±1. The HH entering the final states (12) that con- 
tribute to the result for ^He are listed in Eq. (A. 8). The required ME are sums 
of products of hyperangular ME and hyperradial integrals. The hyperangular 
ME entering the calculation are listed in Eq. (A. 10). The hyperradial integrals 
are of the form (14), they can be calculated analytically and expressed (see 
e.g. [20]) in terms of the hypergeometrical functions. When the model (10), 
(11) is adopted for the bound state, one obtains the following final expression 
for the El strength function: 

di^(El) ^ ^2^^3 (-1)'+^- [alF^iy,)F^iyJ)+a,F,iy,)F,iyJ)] 

Here E is the continuum energy, £"0,1 = (^'«o,i)^/(2m), where Kq, Ki and c are 
defined in Eq. (10). The constant D is 

2 \2m) {A-2)A' 
and the constants ai and a2 are 

1 /315\2 ^ ^ ^ ^ /9OO9V3V, 

ai = - I ^ 1 (l + 3woo + VW00W20), Qi2 =1 I {W20 + W21). 

The coefficients w;ooi ""^20) and W21 entering here are the weights of various HH, 
see Eq. (11). The quantities yo,i are defined as yi = E/{E + Ei), and, finally, 

FM = F ^, 4, 1/) , F2{y) = yF ^, 6, |/) , (16) 

where F{a, f3,'~f, z) is the standard hypergeometrical function [20]. The func- 
tions (16) represent the contributions from final states with K = 1 and K = 3, 
respectively. When the hyperradial function (9) is used the result is obviously 
obtained from Eq. (15) by retaining only the term with i — j — and re- 
placing with 2k,o. At — > the strength function (15) exhibits the typical 
E^ three-body behavior at threshold. (In reality "dineutron" correlations may 
influence the spectrum in the threshold region which is beyond our present 
consideration.) 

In Fig. 3 experimental data for the ^He strength function from Aumann et 

al. [14] (Pb target, i?iab = 240 MeV/nucleon) are compared to our model 
strength functions calculated with Eq. (15). The comparison is done both for 
the energy behavior of the spectra. Fig. 3(a), and in absolute scale. Fig. 3(b). 
From the comparison in Fig. 3(a) one sees that the shape of the spectrum 



9 



is well reproduced by our model, especially in the case of the most complete 
version ^3, see Table 1. A shift of the maximum to lower energies in the latter 
case is related to the inclusion of the contribution from the ground state HH 
with K — 0. Maxima of the contributions from the HH with K — 2 are shifted 
to higher energies since these HH lead in part to excitation oi K — 3 final 
states. 

An interesting feature of the calculated strength functions seen in Fig. 3(a) is 
that when one passes from the hyperradial function (9) to the function (10) 
retaining the same hyperangular dependence (e.g. going from to ^^2) the 
shape of the spectrum changes very little. This is explained by the above men- 
tioned fact that only large p values contribute sizably to the result. For such p 
values only the first, longer range exponential in the hyperradial function (10) 
survives leading to just the same energy dependence as the model (9) (but to 
very different magnitudes of the cross section). Indeed, if one plots a figure 
similar to Fig. 1(b) but retaining only the contribution from the first expo- 
nential to the ME one obtains the integrands extremely close to those shown 
in Fig. 1(b). The same holds true for the K = 3 final state. This means that 
the i — j = term dominates Eq. (15) at not too high energies. The other 
terms are suppressed by the factor [{Eo + E)/(E-i + E)]^^/"^. At higher energies 
the model \&2 leads to a strength function which decreases faster than that for 
the model The internal part of the hyperradial function contributes more 
to the results here, and in the case of the two-parameter model this internal 
part is smaller than for the one-parameter model. 

Comparing the magnitudes of the calculated strength functions with experi- 
ment in Fig. 3(b) one sees that the ground state model ^1 (see Table 1) leads 
to a very low strength function. The results are considerably improved when 
passing to the model \1'2 while the model ^3 leads to a further improvement 
and compares reasonably well with experiment. These results may be com- 
mented as follows. Underestimating the pims value the one-parameter model 
^1 underestimates simultaneously all sizes in the system including the dis- 
tance Tc between the core and the CM, see Table 1. Furthermore, one needs 
to take into account that the El sum rule (5) strictly preserves its value if 
one replaces true final states with a complete set of no interaction final states. 
So, in our case the total strength is determined by the Tc value, and thus it 
is natural that the above model leads to a too low total strength. While the 
two-parameter WF, ^2, reproduces correctly the rms matter radius it still 
gives an Tc value which is somewhat low. For this WF, which includes only 
one HH with = ly, we have the simple relation (x^) = (y^) = (p^)/2. How- 
ever, this relation breaks down, due to off-diagonal terms in the expectation 
values, when more HH are included, as in ^3, and the geometric properties of 
the state change slightly leading to a higher Tc value. 

In accordance with the above discussion, at not too high energy the main 
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effect of passing from tlie one-parameter model (9) to tlie two-parameter 
model (10) consists merely in the multiplication of the strength function by 
c^/{2kq) ~ 4 which can be seen from comparing the strength functions for 
models *i and in Fig- 3(b). This is in principle similar to the case of 
photodisintegration of deutrons [21]. The cross section, calculated with the 
zero-range Yukawa ground state wave function (2ft;o)^/^ exp(— Kr)/r and the 
no interaction final state, is lower than the experimental one. The main effect of 
the finite range correction consists in an increase of the asymptotic constant in 
the ground state WF leading to a reasonably good comparison with experiment 
at moderate energy. 

In Table 2 the integrated strengths obtained are listed along with the exper- 
imental strengths and those obtained from a microscopic three-cluster calcu- 
lation. The strengths are integrated up to = 4 MeV, 9 MeV, and infinity. 
Note that for the first of these energies the ratio of the result obtained with 
to that obtained with is about c^/(2ko) as it should be. Our total strength 
proved to be very close to that obtained using the microscopic bound state 
WF. These results on the total strength provide us also with a check of the 
strength function calculation. Integrating the strength function up to infinite 
energy one should reproduce the value of the sum rule (5). 

Finally, we calculate the EMD excitation energy spectrum using Eq. (1). The 
virtual photon spectrum was obtained for a ^He beam with energy 240 MeV/A 
striking a Pb target. In Fig. 4 the resulting distribution is compared with the 
excitation energy spectrum measured by Aumann et al. [14]. In Table 3 the 
total ^He cross sections obtained with our model are compared with experi- 
mental estimates for two different beam energies, 37 MeV/A and 240 MeV/A, 
measured at NSCL [22] and GSI [14] respectively. The theoretical cross section 
has been calculated by integrating the spectrum up to the energy £■ = 10 MeV, 
corresponding approximately to the maximum energy measured in the exper- 
iments. The observed difference should be a direct measure of the importance 
of FSI since the virtual photon spectrum peaks at low energies and this is 
the region were FSI is supposed to play a role. We find that our model gives 
~ 35% lower cross section than experiment for both energies. 



4 Conclusion 

A three-body model describing EMD of two-neutron halo nuclei, which allows 
studying the El strength functions without precise knowledge about initial and 
final states, has been developed. The model leads to an analytical expression 
for the strength functions. Our ground state wave functions reproduce both 
the true three-body asymptotics (determined by the binding energy) and the 
correct size of the system. A complete set of three-particle no interaction final 
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states is used. We found that the large distance asymptotics of the ground state 
determines the shape of the strength function, while the size of the ground 
state governs the asymptotic constant and thus the magnitude of the strength 
function. We pointed out that in some cases the shapes of the El strength 
functions for different two-neutron halo nuclei are approximately related to 
each other. 

We have used ^He as a test case and made an extensive comparison with 
experimental data [14,22]. We found a good agreement concerning the shape 
and peak position of the El strength function and a reasonable agreement 
concerning its magnitude. 

A remarkable feature of halo nuclei is that the El strength is concentrated at 
low energies even without any low-lying resonant state. This peculiarity is due 
to the low binding energy and large size of the initial state which is demon- 
strated by the present calculation with no final state interaction included in 
the model. 

In a forthcoming paper, see also [23], we shall apply the present approach 
to investigate the El strength functions of Borromean ^^Li and ^^Be nuclei. 
Common for these nuclei are large uncertainties concerning their microscopic 
structure. In the ^^Li case the El strength function is not so well known 
experimentally as in the ^He case, and for ^^Be the EMD process has been 
measured only very recently [24]. 

As to the FSI effects, we noted that the total El strength will not change when 
replacing a complete set of no interaction final states with such a set of true 
final states. Therefore the only effect of including FSI will be a redistrubution 
of the strength leading to a somewhat higher strength at low energy. Further 
studies on this point would be of interest. 
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A Appendix 



A.l Coordinate sets and hyperspherical harmonics 



We need to perform the calculation in the CM subsystem. The adopted Jacobi 
coordinates are 



y = 



l2{A-2) fn+f2 



(A.l) 



where ri, r2, and fc are positions of the valence nucleons and the core. We use 
the related hyperspherical coordinates {pjfl^}. The hyperradius p determines 
the size of a three-body state: 



P 



+ = {fi - Rcu^f + (r2 - ^cm)' + {A- 2) (re - ^cm)'. (A.2) 



The five angles {Vl^} include usual angles {9x,(f>x), {(^y,(f>y) parametrizing the 
unit vectors x, y and the hyperangle 9 defined by the equalities 



X — psinO, y = pcos9, 
where < 6 < 7r/2. The volume element dfl^ entering (7) is 

dQ5 = sin^^cos^^d^dxdy. 
The HH have the explicit form 



r^W(^5) = #(^)ili(x,y), 



where 



>li(x,y) 



(A.3) 



(A.4) 



(A.5) 



(A.6) 



and Yim are spherical harmonics. In Eq. (A.5) the hypcrangular functions are 

^^'''(^) = N^-^" sin'- ^cos'^ eP^-+^/^'^y+^/^{cos2e), 

where P^'^{x) are the Jacobi polynomials (see, e.g., [20]), n = [K — l^ — ly)/2, 
and N'p are normalization constants. 



K 



2{n\){K + 2){n + U + ly + l)\ 



\ V{n + + 3/2)r(n + ly + 3/2) ' 
The HH (A.5) are orthonormalized using the volume element (A.4). 
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A. 2 El transition matrix elements 



The explicit expressions for the HH entering our ground state are the following 
(see (A. 3), (A. 6) for notation): 

'IT 



= — cos2^roo°(x,y), (A.7) 



The explicit expressions for the final state HH contributing to our strength 
function are the following: 



r?U(fi5) = -^cosey°^jx,y), 



ri?M,W = -|^siney/^jx,y), 



niM, (^^5) = -^(4 COS 29 + 1) cos e (x, y), (A.8) 
nUA^^) = ^(4cos2^ + l)sin^yii^Jx,y), 



32 

^IIm,{^5) ^ ^cos^ 9sm9Yll,^{±,y), L = 1 and 2. 

V OTT 



The El operator, Eq (4), can be written as 



M{El,f,) = -eZ,^ ^J_^^ pcos9Y,^{y). 

The ME 

{Kf lyf Lf S, Jf = 1, M/ = /i| cos^yi^(y)|Xi k hi U S, Ji = 0) (A.9) 
of the hyperangular part of this operator, between basis states of the form 



JM 



that include the HH (A.7) and (A.8) are required. We denote these ME as 

{Kf Ix lyf Lf S\0\Ki Ix lyi Li S). 
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The non-zero ME are 



(10110|O|00000) = 

(30110|O|20000) = 
(31011|O|21111) = 
(31221|0|21111) 



1 



2V27r' 



4V27r' 
1 



1 



(10110|C>|20000) 

(11011|O|21111) 
(31211|C>|21111) 



1 



4^27?' 
1 

1 



(A.IO) 



4V37r 



4A' 



The calculation with respect to spin variables and angles x, y is done using 
the formula (see, e.g., [25]) expressing the ME (A. 9) in terms of reduced ME in 
the X, y subspace and then in terms of reduced ME in the y subspace. Finally, 
note that the results obtained in this section are useful not only for ^He but 
also when studying other two-neutron halo nuclei. 
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Fig. 1. (a) The model hyperradial function x^'^\ Eq. (10) (solid line), compared 
with the three predominant hyperradial functions in expansion (6) obtained from a 
microscopic three-body calculation of ®He [17]. All functions have been normalized 
to unity, (b) Typical integrands entering the El transition ME (14). This plot 
illustrates the important range of p values. 

1.4r 




Fig. 2. Comparison of scaled, experimental strength functions for He and Li. 
The experimental data are: (circles) ^He - Aumann et al. [14], (triangles) ^^Li 
- Sackett et al. [26], (squares) ^^Li - Shimoura et al. [27]. The dotted curve is 
the two-body strength function from Ref [10,11]. The solid line is the three-body 
strength function from Eq. (15). 
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E (MeV) E (MeV) 

(a) (b) 

Fig. 3. A comparison of our strength function with experimental data for ^He, Au- 
mann et al. [14], both in arbitrary units (a) and in absolute scale (b). The analytical 
strength functions are obtained with WFs (dotted), "if 2 (dashed) and ^'3 (solid), 
see Table 1. 

150| . ■ 1 




E (MeV) 

Fig. 4. The ^He EMD excitation energy spectrum obtained at 240 MeV/ A with 
a lead target. The data points are extracted from Ref. [14] while the solid line is 
our result obtained by multiplying the strength function (15) calculated for the ^'3 
model with the spectrum of virtual photons. The dashed line is the same curve 
corrected for the LAND detector response and efficiency. 

Table 1 

The model WFs for ^He. The parameter kq reproduces the two-neutron separation 
energy. The last three columns give the hyperradius Prmsj the matter radius i^rms 
and the distance between the core and the CM r'c. 

WF Kq Ki Wqq W20 W21 Prms Rims 

(fm-1) (fm-i) (fm) (fm) (fm) 

*i 0.2163 — 1 3.27 1.79 0.67 

*2 0.2163 0.5370 1 5.37 2.50 1.10 

*3 0.2163 0.5370 0.05 0.80 0.15 5.37 2.50 1.20 
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Table 2 

The ^He integrated strength dE' dB (El) / dE' , where E is the continuum en- 
ergy, for different models and experimental data. The last column corresponds to 
E = oo and equals to the result of Eq. (5). All values are given in e^fm^. 





E = A MeV 


E = 9 MeV 


Total B(E1) 


Aumann et al [14] 


0.59(12) 


1.2(2) 






0.079 


0.18 


0.43 




0.29 


0.63 


1.14 


^3 


0.42 


0.83 


1.38 


Realistic WF|3 


0.71 


1.02 


1.37 



Danilin et al. [9] 



Table 3 

The ^He EMD cross section obtained with a lead target at different energies. The 
theoretical values are obtained by integrating (1) up to 10 MeV. 



Energy (MeV/ A) 


Exp. 




This work 




CTEMD (mb) 


Ref. 


CJEMD (mb) 


37 


830[] 


[22] 


551 


240 


520 ±110Q 


[14] 


333 



^ 80% of the incident ^He ions were in the range 28-52 MeV/A. A 
Glauber-type calculation estimated that 60% of the total inelastic 
cross section were due to EMD. 

^ The EMD cross section was obtained by subtracting eikonal 
model cross section [28] (for nuclear excitations) from the mea- 
sured inelastic excitation cross sections. 
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